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ABSTRACT 


A CO? laser was used to generate reduced-density channels in various gas mixtures 
of ammonia and nitrogen. Interferometers were used to record the changing density re- 
sulting from the NH3 absorption of the CO2 laser radiation. One method used to de- 
termine the characteristics of the channel was based on the assumption that the resulting 
density profile was Gaussian shaped. The second method used the Abel [ntegral Trans- 
formation, requiring no pre-conditions on the density profile except that it was 
cylindrically symmetric. While used extensively in plasma spectroscopy, this technique 
is not generally well known for analyzing interferometric data. The results show that for 
a fixed laser energy long, shallow reduced-density channels were formed in gas mixtures 
of low ammonia concentration and short, deep channels Were formed in gas mixtures 
of high ammonia concentration. These results qualitatively agree with an earlier exper- 
iment in which gas mixtures of nitrogen and sulfur-hexaflouride were used. Both exper- 
iments support the concept of reduced-density channel formation and lay the foundation 
for future studies of relativistic electron beam propagation for application in the Strate- 


gic Defense Initiative. 
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I. INTRODUCTION 


The Strategic Defense Initiative (SDI) was conceived with the goal of making nu- 
clear weapons obsolete. One of the many areas being investigated in the pursuit of thus 
goal is the use of intense high energy electron beams [Ref.1]. In the investigation of 
these beams, it was found that interaction with the atmosphere reduced electron beam 
propagation. To improve beam propagation, a method was sought to reduce the at- 
mospheric density along the path of the beam. This method, Reduced-Density Channel 
formation, consisted of using a laser to pre-condition the atmosphere bv creating a re- 
gion of reduced density along the intended path of the electron beam. 

Studies of Reduced Density channels were conducted prior to the establishment of 
fmemotratesic Wefense Initiative [Refs, 2, 3]. These studies led researchers at the Los 
Alamos National Laboratory to conduct experiments in channel formation. Two groups 
of experiments were conducted in which variable laser energies of up to 100 joules were 
incident on target cells containing gas nuxtures. One group of experiments used target 
cells containing mixtures of nitrogen and sulfur-hexaflouride [Refs. 4 through 7]. The 
second group of experiments which used target cells containing various nuxtures of mi- 
trogen and ammonia are discussed here. A schematic of both experiments 1s shown in 
Figure |. 

A COQ2 laser beam incident on a target cell containing a nuxture of NH3 and N2 
creates a region of reduced density when the NH3 partially absorbs the laser enegv. The 
Characteristics of the reduced-density channel are determined by analyzing 
interferograms taken after the laser was fired. 

The CO2 laser used in this experiment was tunable and produced laser pulses up to 
300 joules. The laser pulses exhibited a time profile shown in Figure 2 and had a 
guassian spatial profile. Approximately 33 percent of the beam energy passed through 
the aperture and energized the gas in the cell: that is a maxunum of 100 joules was 
available for channel formation. As the beam traversed a cylinder filled with a gas 
mixture of NH3 and N2, a region of reduced densitv along the longitudinal axis was 
formed. Three interferometers were positioned at different locations along the cvlinder 


meee Ficure |), 


Pictures of the interference patterns, interferograms, were taken to record the eileen 


the laser beam on the gas. These pictures were then digitized and the infornnationesises 
to determine the depth and width of the resulting channel. Two techniques were used 
to determine the charactensties of the densi -cuammer 

The first technique, the Gaussian method, assumed that the density profile resulting 
from laser absorption had a Gaussian shape. The second technique developed to deter- 
mine the channel characteristics required the assumption that the density profile was 
cylindrically symmetric. The Abel integral transformation, a well known plasma diag- 
nostic technique, was adapted for use with interferometric data. This adaptation is be- 
lieved to be one of the first applications of the Abel transformation to interferometric 
data. 

This thesis details the procedures used to extract reduced-density channel charac- 
teristics from interference patterns recorded during the experiment and constitutes a 
thorough analysis of all the data from the second group of the Los Alamos experiments. 

These results, described in Chapter 6, will indicate: 


e Fora fixed energy level, increasing the concentration of NH3 results in deeper and 
shorter reduced-densitv channels. 


e Fora Gaussian shaped interference pattern, the Gaussian method of analvsis pro- 
vides a quick and accurate means of determining the reduced-densitv channel 
characteristics. 

These results are supported by the N2-SF6 experimental results and will be used to studv 


the propagation of intense electron beams through the atmosphere. 


H. DESCRIPTION OF THE REDUCED-DENSITY CHANNEL 
EXPERIMENT 


A. EQUIPMENT 
A schematic of the experiment is shown in Figure |. A tunable electron beam con- 
trolled CO2 laser with an 11x11x100 cm? active volume produced pulses of 10.6 um ra- 
diation at energies of up to 300 joules having a 3.3 ws pulse width. Upon leaving the 
laser, the radiation experienced multiple reflections as 1t passed through a nitrogen filled 
cylindrical tube 10 m long. The reflections served to increase the optical path, hence the 
laser beam traveled approximately 35 m before exiting the tube and entering the test cell. 
The increased length aided in focusing the beam; the nitrogen in the cylinder prevented 
the beam from being scattered bv particulate matter. Since nitrogen does not absorb 10.6 
um radiation the laser beam approached the apperture without anv significant degrada- 
tion. One-third of the initial laser energy passed through the aperture and was incident 
upon the gas at the entrance of the target cell. A typical 60 joule pulse entering the target 
cell had an estimated peak on-axis fluence of 50 Joules per c?. Entering the cell with 
a full width at half maximum of ~10 mm, the beam was aligned so that the peak of its 
Gaussian profile was along the longitudinal axis of the cylindrical test cell. [Ref. 8: pp. 
4-5}. 
The test cell, a cylinder 36 cm in diameter and 250 cm long, contained a mixture of 

NH3 and N2. Experiments in the test cell, were conducted on the following mixtures: 

Se vorr iNI13 and 595 torr N2 

e 15 torr NH3 and 585 torr N2 

e 50 torr NEH3 and 550 torr N2 

e 200 torr NH3 and 400torr N2 


The test cell was evacuated before each experiment. Subsequently it was filled to the 
desired pressure of NH3 and then brought to atmospheric pressure (600 torr for Los 
Alainos) with pressurized N2. 

The NH3 molecules absorb 10.6 um radiation [Refs. 9, 10 and I]]. As the beam il- 
luminated the cell, the absorbed energy was transferred to the other molecules of the gas 
mixture through collisions. The adiabatic expansion of these excited molecules produced 


aeeduced density chantiel which was centered on the longitudinal axis of the test cell. 


An tn depth description of the adiabatic expansion process is developed in reference 
8 pages 5 through 10. Included in that reference are calculations which determine the 
required absorbed fluence to produce a unform half density channel 2 min length. While 
the values obtained therein apply to 600 torr N2. they apply to the mixtures used in this 
experiment duc to the prominence of the N2 contained in each gas mixture? Ji/shoruume 
fluences provided by the CO2 laser were sufficient to produce 2 m long channels which 
were | cm wide and had 0.1 times the ambient densitv. 

In order to observe this reduction in density, windows were installed along the walls 
of the cell. These windows permitted the use of optical interferometers, one of whose 
beamis traversed the test cell. The density profile of the gas nuxture at different locations 


was then determined by analvzing the interference patterns. 


B. INTERFEROMETERS 

Three interferometers were installed at different positions along the test cell (See 
Figure 1). These Mach-Zehnder type of interferometers spht the He-Ne laser (633nim) 
into two separate beams. One beam traversed across the path of the CO2 laser beam 
and the other beam traversed around it. Figure 3 1s a photograph of the pattern resulting 
from interference of both beams prior to the energizing of the CO2 laser and subsequent 
channel development. This pattern is due to non-alignment of the murrors used in the 
interferometer [Ref. 12]. Filling the test cell with the NH3 and N2 mixture and enersizme 
the CO2 laser developed a phase shift between the two beams. This phase shift was due 
to a decrease in the optical path length of the beam passing across the path of the CO2 
laser. Figure 4 is a photograph of the resulting pattern. 

Photographs were taken of each interference pattern observed. These photographs, 
interferogranis, Were processed in order to deternune the density profile. The processing 
required photographic enlargment, digitizing and numerical analysis. The following sec- 


tion 1s devoted to describing terms used in analyzing the interferograims. 


C. INTERFEROGRAMS 

An interferogram is a photographic record of a pattern produced by an 
interferometer. In this experiment three separate Mach-Zehnder type interferometers 
were used at different locations to determine the density profile of a test cell containing 
a gascous mixture of NI13 and N2. The following characteristics are used to describe 


the pattern in an interferogram. Figure 5 is a sketch with these terms clearly annotated: 


Top and Bottom of interferogram. 

The edge of the interferogram towards which the peaks of the fringe lne are 
displaced is the top. The bottom of the interferogram is the edge directly opposite 
The top. 


Pc felerence mark. 

In order to establish a scale between the interferograms, a 2 centimeter refer- 
ence mark was devised and located in the field of view of the camera. This device, 
a periodic 2 em square-tooth pattern was used for each interferogram. Portions of 
Maerdevice are visible im migures 3 and 4. Im Figure 5 tlie device is shown at the top 
of the sketch. The distance from the right-most corner of one tooth, to the right- 
most corner of the next tooth was 2 cm. 


Prince liye: 

The interference pattern showns in Figures 3 and 4, also sketched in Figure 5, 
consist of alternating regions of black and white. Although either region could be 
used to define a fringe line, the black region 1s more clearly defined visually. Thus 
throughout this description, a fringe line refers to a line in the black region which 
is halfway between two white regions. Figure 5 shows a sketch of two typical fringe 
lines. 


Fringe number 

The fringe number was arbitrarily chosen to be the number of fringe lines from 
the bottom of the interferogram. In Figure 4, the third fringe 1s found bv starting 
at the first black region at the bottom of the interferogram. The fringe line in this 
region is fringe number one. The next black region encountered as one moves from 
bottom to top 1s fringe number two. contains fringe number two. The next region 
contains fringe number three and so on. 


Zero fringe line. 

A local term used to denote a straight line passing from the left-most edge of 
a fringe line to the right-most end of the same fringe line. This line corresponds to 
the interference pattern which would be present if there were no fringe displacement 
(1.e., no absorbing gas in the test cell). Two of these zero fringe lines are annotated 
Merigure 5, 


Fringe displacement. 

The distance measured in fringe lines above the zero fringe line. It is due to 
the decreased optical path length created when the gas mixture absorbs 10.6 sm 
radiation. To determine fringe displacement, one must find the distance between 
successive zero fringe lines. A ratio of these numbers yields a scale factor from 
Which any distance measured along the vertical can be converted to a fringe value. 
In Figure 5, fringe number one has Zero fringe displacement at the far left and the 
far right. The peak of this fringe line is approximately 1.25 fringe lines. 


Vertical reference line. 

‘\ Vertical reference line was drawn connecting the peaks of the fringe lines. 
This line was paralle! to the vertical faces of the square tooth but was not neccs- 
sarily perpendicular to the zero fringe line. 


D. SUMMARY OF EXPERIMENT 

In this chapter the experiment and its recorded results are described. The next step 
in trving to determine the densitv profile and subsequently the characteristics of the 
reduced-densitv channel, 1s to interpret the interferograms. This interpretation is con- 


ducted in the next chapter. 


Wi. DATA REDUCTION 


Photographs of the interference patterns were taken shortly after the CO2 Laser 
energized each gas mixture. These interferograms provided a record from which the 
fringe data, the relationship between fringe displacement and displacement along the 
zero fringe line, would be obtained and the on-axis densitv determined. Once the fringe 
data was obtained the density profile and subsequently the reduced-densitv channel 
characteristics could be determined. This chapter describes how the fringe data was ex- 
tracted from the interferogram and placed into data files. These files will be analyzed 
later by the Gaussian and Abel techniques to determine the characteristics of the 


reduced-density channel formed in each gas mixture. 


A. PHOTOGRAPHIC ENLARGEMENT 

Determining the fringe data required making measurements directly on the 
interferogram. To facilitate this procedure, each original interferogram, a photograph 
measuring approximately 3 cm by 3 cm, was enlarged. The enlargement process was 
accomplished using a specially modified 35 mm camera fitted with a zoom lens. The 
mage observed in the field of view of the camera was converted to an electronic signal 
and simultaneously passed to a video screen and Tektronix photo-copier. While moni- 
toring the video screen, the focus of the camera was adjusted until a clear image of the 
2 cm reference marks and interference pattern filled the video screen. The ‘copy’ button 
of the photo-copicr was then depressed, resulting in an enlarged photograph of the 
interferogram. Some interferograms were not capable of being enlarged with the 2 cm 
reference marks visible. For those intcerferograms, a scribe was used to mark the refer- 
ence points on the orignal photograph. The scribe marks enhanced the image of the 
reference marks and made them visible in the enlarged interferogram. 

Each enlarged photograph therefore contained a set of marks which measured 2 cm 
in the original interferogram. The distance between these marks varied with each 
interferogram as a result of the focusing adjustments made during the enlargement 
process. Thus, all distances could be scaled to the comunon 2 cm reference mark and 


continuity was maintained between photographs. 


Having enlarged each interferogram. the next step in determining the density profile 


and subsequently the characteristics of the reduced-densitv channel was to digitize the 


enlarged photographs. Before digitizing the interferograms, a one-time check of the 
combined effects of digitizing and enlarging was conducted by comparing an original 
drawing to an enlarged copy of the same drawimg. A Ine was drawn on a piece of graph 
paper and eight easilv identified points were then digitized (See Figure 6). The drawing 
was then enlarged, using the process outlined above, and the same eight points on the 
enlargement were digitized. The orginal and resulting coordinates are listed in Table 1. 
A review of this Table illustrates that the maximum difference between the original and 
the enlarged copy are within .04 cm of each other. The mean difference between the or- 
iginal and the enlarged copy was no greater than .02 cm for both the vertical and honn- 
zontal components. It will be shown in chapter 6 that the errors generated bv the 
enlargement process did not significantly effect the results. The next step in determining 


the density profile was to prepare the interferograms for digitization. 


B. DIGITIZING PREPARATIONS 

In preparation for digitization, each enlarged interferogram was reviewed and the 
following items identified by marking on the photograph. Figure 7 1s a sketch of an 
inteferogram prepared for digitizing. 


e 2cm. reference marks. 
The 2 cm reference mark was identified by drawing a line from the lower right 
corner of the left-most square tooth, Rl, to the same corner of the next square 
tooth, R2. | 


* lyrince lines to be dicitized. 
Two fringe lines per interferogram were chosen to be digitized. These fringe 
lines Were normally picked based on their clarity and proximity to the center of the 
interferogram. Three fringe lines labelled Fl, F2, and F3 are shown im Figures 


© Verticalirelcrence dine: 
A vertical reference line was drawn connecting the peaks of the fringe lines. 
This line was parallel to the vertical faces of the square tooth but was not neces- 
sarily perpendicular to the zero fringe line. | 


* Zero innge lines: 
Zero fringe lines were drawn for many of the fringe lines of each interferogram. 
These lines provided the basis for determining the fringe displacement and subse- 
quently the on-axis density. Three zcro fringe lines have been drawn in Figure 7. 
Having prepared each interferogram for digitizing, the next step was translating the 
photographie information available from the interferogram to a form which would allow 


humerical processing. 


C. DIGITIZATION 

Converting the photographic image of the interference pattern into a form which 
allows mathematical computation 1s the process of digitizing. 

The first step in the process was to establish a rectangular coordinate system for the 
interferogram. This was accomplished by designating the left-most point on the fringe 
line as the origin. This point is labelled P(1) in Figure 7. The vertical axis was made 
parallel to the vertical reference line of the interferogram and the horizontal axis was 
chosen perpendicular to this. With this reference system established any point on the 
interferogram could be designated by a vertical and horizontal component. The 
HP-9111A Graphics Tablet was used to assign numerical values to these components. 
The numerical values are in digitizing units, where | digitizing unit 1s equal to 0.025 
millimeter {Ref. 13]. The Graphics Tablet can be thought of as an electonic grid with its 
own rectangular reference system and units. The wand used with the Tablet had a spring 
loaded tip. With the tip in contact with the Tablet and the spring compressed, the posi- 
tion of the wand relative to the reference system of the Graphics Tablet was registered. 
This information was then passed to the HP-85 computer and stored in a data file. 

Securely fastened to the Tablet using adhesive tape, the interferogram was ready to 
be processed. This consisted of digitizing the following characteristics of the 


interferogram. 


wa 2ici reference mark. | 
Digitizing the right-most corner of one square tooth and then the same corner 
of the next tooth allowed the distance, in digitizing units, to be deternuned. This 
distance corresponded to 2 cm on the interferogram. Refer to points labelled Rl 
and 2 on Figure 7, 


e Distance between zero fringe lines. 

Pigeon tie Mitersection of the vertical reference line and a zero fringe lime 
and the vertical reference line with different zero fringe line allowed a scale factor 
to be determined. The distance between these two digitized points represented the 
number of digitizing units per fringe. Refer to points R3. R4 and R5 of Figure 7. 
Having previously determined the number of digitizing units per 2 cm, a value for 
the number of fringes per cm was determined. 


e Fringe line. 

All fringe lines, two per interferogram, Were digitized in the same manner. 
Starting at the left-most point of the fringe line, point P(1) Figure 7, and proceeding 
towards the right. approximately 35 points per line were digitized. (Sce fringe line 
Fl of Figure 7.) The last point digitized was the right-most point of the fringe line, 
P(N1). 


As the points on the interferogram were digitized thev were being stored in a data 
file. Each digitized fringe line, along with the digitized reference points. was recorded in 
a separate file. 

Each file contained enough information to reproduce the fringe line as it had been 
digitized. In order to regain the relationship between fringe displacement and distance 
along the zero fringe line, some data manipulation was needed. This manipulation 
consisted of correcting the digitized fringe line for the angular difference between the 
zero fringe line and the horizontal reference line. At this point it is important to recall 
that the fringe value is determined by the distance the fringe line is away from the zero 
fringe line-not the horizontal reference line. The following outlines the method used in 


Writing a computer program which converted digitized data into fringe data. 


D. FRINGE DATA 

The fringe displacement as a function of distance along the Zero fringe line defines 
a rectangular coordinate system. In order to obtain this relationship, a coordinate 
transformation was made. This transformation required aligning the horizontal axis of 
the digitized coordinate system with the zero fringe line. The process outlined below 
follows closely the development of a computer program which performed the transfor- 


mation. To assist in this development, the following terms are defined. 


e J is anumber assigned to the point digitized. J=1 for the first point digitized, J=2 
for the second, and so on. 


¢ F is the number of zero fringe lines between the fringe lines digitized for ima 
distance calibration. In Figure 7, the number ts | if R3 and R4 were digitized. F 
equals 2 1f R3 and R5 were digitized. 


e V6 is the angle, in radians, between the zero fringe line and the horizontal axis. 
¢ D1 is the distance, in digitizing units, between the 2 cm reference marks. 


e D2 is the distance, in digitizing units, between zero fringe lines digitized for a ref- 
erence. 


e Nl is the last point of the fringe line digitized: the right-most point. 
¢ XI is the horizontal component of a point in digitizing units. 
e Yl is the vertical component of a point in digitizing units. 


e A(J,1) 1s the horizontal component of a digitized point J. 


A(J,2) 1s the vertical component of a digitized point J. 


In Figure 8, one of the many fringe lines of an interferogram has been sketched and 


annotated to show the g 


o— 
ter 


eometry of the coordinate transformation. Having identified the 


geometric relationships in the interferogram, the next step in generating fringe data was 
referencing all digitized points to the origin. This referencing was accomplished by 


subtracting the coordinates of the origin from each point digitized. As shown: 
X1 = A(J,1) — A(1,)) (oe) 
Y1 = A(J,2) — AC,2) (ge2) 


The next step was to determune the angle between the Zero fringe line and the hori- 
zontal axis. This was accomplished by using the first and last points of the digitized 
fringe line. The ratio of the difference of vertical components to horizontal coniponents 


was the tangent of the desired angle, V6. 


X¥6 = A(N1,1) — A(1,1) (3.3) 

owe 12) <wi(1,2) G4) 
a7 6 : 

V6 = ARCTAN > (3.5) 


With this angle determined, the next step was to relate the horizontal value to a 
value along the zero fringe line. Since the horizontal component was the projection of 
the zero fringe line onto the horizontal axis, the distance along the Zero [ringe line, R$, 
was found by: 

AG 


§ EN (3.6) 


Converting this into centimeters vields the horizontal component of the fringe data 


R8&(J) = R8x2x D1 (3.7) 


The final step in determining the fringe data was to find the fringe value ofa vertical 
component. This was done bv finding the contribution to the vertical component made 
by the zero fringe line and subtracting it from the digitized vertical component. 


The vertical contribution of the zero fringe line ,Y8, was determined from: 
¥8 = V1 tan( V6) (aio) 


This value was subtracted from the vertical component of a digitized point and 


finding its projection along the horizontal axis. 


D3 =(Y1— ¥8) cos(V6) (3.9) 


The fringe value was finally obtained by multiplying D3 by the ratio of fringes to 


digitizing units. 


Pe Fx D3 ‘ 
YS(/) = ae (3.10) 
Hence the fringe data originally observed on the interferogram was converted into 


a numeric form, R8(J) , Y(J) which allowed mathematical processing. 


E. SUMMARY 

Digitizing the interferogram was necessary to allow mathematical processing by the 
computer. In the process, determining the relationship between the fringe displacement 
and the zero fringe line was required. Having accomplished this, determining the density 
profile and subsequently the characteristics of the reduced-density channel, using the 
computer was possible. This determination was conducted using two separate tech- 


niques which are addressed in the following chapter. 


IV. GAUSSIAN ANALYSIS 


A Reduced-Density Channel was generated by firing a CO2 laser beam into a test 
cell filled with a gas mixture of N2 and NH3. Mach-Zehnder tvpe interferometers were 
positioned at different locations along the test cell and the resulting interference patterns 
were recorded on photographic film. These records, interferograms, Were analvzed using 
Gaussian Analysis and Abel Inversion techniques. This chapter 1s devoted to the devel- 
opment of the Gaussian technique and the numerical manipulation required to obtain 
the characteristics of the reduced-density channel from the digitized data files using the 


Gaussian technique. The Abel technique is addressed in the next chapter. 


A. THEORY 
The Gaussian technique used to determine the characteristics of a reduced density 
channel was originally developed by Czuchlewski for his experiment with a test cell 
containing N2 slightly doped with SF6 [Ref: 8: pp. 11-14]. This technique will be pre- 
sented by first considering a uniform square channel and obtaining a relationship be- 
tween the on-axis density, fringe displacement and full width at half maximum. The 
square channel case will then be generalized to cover a non-uniform channel. The result 
will be a relationship similar to the square channel case, but applicable to the 
interferograms obtained in the reduced-density channel experiments. In this develop- 
ment a link between the Gaussian and Abel techniques will be identified along with key 
values used in both techniques. These key values are listed in Table 2. Having deter- 
mined a relationship between the on-axis density, fringe displacement and full width at 
half minimum for a non-uniform channel, its application to the stored data files will be 
described. The results of the Gaussian Analysis for each fringe line are presented in 
Tables 3 through 6. 
1. Square-Channel 
Consider a hypothetical situation in which one leg of the interferometer is 
passing through a square channel. This channel and the resulting interference pattern 
are shown in Figure 9, Designating the width of the channel as I. the optical path dif- 
ference, OPD, between the arms of the interferometer can be related to the index of re- 


fraction in the channel and in the reference beam by: 


OPD=(n.—1,) xT (4.1) 


Wanere 
* 7. is the index of refraction im tie chammae 


e nv, is the index of refraction in the reference arm of the intercromercr. 


One fringe shift equates to a phase change of 2x. The fringe shift identified in 
Figure 9, AF. is the optical path difference divided by the wavelength, 4. Substituting this 


in equation 4.1 vields a relationship between fringe shift, AF, and channel width, I. 
AF = (n~ %) (4.2) 


The index of refraction is related to the mean polarizability, a, of a gas in CGS 


notation through the Lorenz-Lorentz relation [Ref. 14: p. 87] 





3 cere 3 Ee-!l (4.3) 


where 
® ¢ 1s the dielectric constant. 
® nis the index of refraction. 
¢ Nis the number of molecules per unit volume. 


® ¢=}I- 


Making use of the approximation that the index of refraction of a gas is nearlv 


unity, equation 4.3 can be reduced to: 


no — | 
tS anN (4.4) 


Rearranging equation 4.4 and solving for the index of refraction yields: 
tZ, 
n= (1 +4204)” (a) 


Applying the binomial expansion theorem to equation 4.5 and discarding the higher or- 


der tems wields: 


H=1+ 2raN (4.0) 





In order to obtain a relationship between the index of refraction and the density, 
a connection between the mean polarizability, « , and densitv must be inade. The perti- 


nent connection is made through the molar refractivitv, A [Ref. 14: p. 88]. 


2 
4 _ RT. {nr —-l) Bes 
=—x ze\. = XX +i 
A 7 X taNy eas (4.7) 
where 


e Ris the ideal gas constant 
e T is the temperature in degrees Kelvin 
e P is the pressure of the gas 


e N,is Avagadros number. 


Solving equation 4.7 for the index of refraction, n, and using only the first order 


terms of the binomial expansion theorem vields: 


Belt a (4.8) 
Recognizing that the density of a gas, p , can be expressed as: 

p= (4.9) 

and solving equation 4.8 for the index of refraction in terms of density yields: 

fe = =e (4.10) 


Ne 


From equations 4.3 and 4.7 one observes that the mean polarizability, «, and the molar 
refractivity, A, are constants. This is confirmed by Born and Wolf [Ref I4: p. 88] as they 
Seamonstrate that the molar refractivity of air varies less than 2 percent between pres- 
sures of 1 and 42 atmospheres. In their development of the Gaussian analysis, 
ee ichiewski ct al., [Ref. $: p .}1], regard the index of refraction as a product of the 


density. p, and a constant. Flence equation 4.10 can be written as: 
n—1l=Kp (4.11) 


where 





a (4.12) 


Substituting equation 4.11 into equation 4.2 vields a relationship between the 


fringe shift and the density of the gas nuxture. 


= (K.0¢ iat eal 
A 


ue (4.13) 
where 


@ p,is the density of the gas in the channel 


e pp, is the density of the gas at ambient conditions 


Since it has been determined that A. and subsequently Kk, is constant, equation 


4.13 can be written as: 


KI (p, ay) 
A 


AF = (4.14) 


Defining 6 as the ratio of the on-axis density to the density at the channel wall, 


allows equation 4.14 to be written as: 


(Oo -1)KTp, 
} 


Be (4.15) 
Hence equation 4.15 provides a relationship which allows one to determine the 
ratio of densities in a square channel test cell by measuring the fringe shift and the full 


width at half maximum. 


AF? 


j=l ae 


(4.16) 
Having developed equation 4.16 for the case of a square channel, the next step 
is to generalize the procedure for a channel which 1s indicative of the interferograms. 
2. Non-Uniform Channel 
Consider the interferogram sketched in Figure 10, This represents a non-uniform 
channel in which the channel density and index of refraction are functions of positians 
(X,Y). The optical path difference at a position along the y axis can be obtamied bilan 


NILE Sralenlome OG wants. 
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L 
OPD(¥) = moss — n JAX (4.17) 


where L is large enough to extend bevond the limits of the channel. Using the definition 
of fringe shift, A F, and replacing the indices of refraction by there equivalence in terms 


of density vields: 





AF(Y) = 


Kp, re 
: [6% - max (4.18) 


Therefore the square channel case, when expanded to cover a general non- 
uniform channel vields an integral relationship between the fringe shift and the ratio of 
densities. 

The method chosen by Czuchlewski to solve the integral [Ref. 8: p. 13] was to 
pick a density profile which was cylindrically symmetrical. The function chosen, 6 (r), 


was a Gaussian function and its characterisitics were similar to the actual 


interferograms. 

2 
§(r) = 1 — (1 — 5,) exp( al (4.19) 
where 


e 6, is the ratio of on-axis density to ambient density. 


e Tis the full width of the channel at the half minimum density. 


In order to solve equation 4.18 the chosen density profile was converted from 
cvlindrical to cartesian coordinates. This conversion allowed equation 4.19 to be written 
as: 


. Soe 
O(X,Y)=1— (1 —0,) expf aie 


] (4.194) 
Upon substituting equation 4.19a for 6 (X,Y) in equation 4.18, one obtains the 


following expression for the fringe shift: 


ee) — —— (1 — d,) exp 


; ee 





Kp, era yf" MBN, 


The integral of equation 4.20 1s Gaussian and its solution allows the fringe shift 
to be obtained from: 


og) Sa | 
B00 = aa] (ES 0) San (4.21) 


Evaluating equation 4.21 at Y=0 vields an expression for the non-uniform 


channel density which closely resembles the relation obtained for the square channel. 
Kp,T 
AF(O) == [5 1*(1- 4.) (4.22) 


Since AF and T are measurable quantities on the interferogram, the ratio of 
on-axis density to ambient densitv can be determined from the following: 


a AF{O) , en wee 
ogee LX, | = taal (4.23) 


Hereafter the 6, will be referred to as channel density. 


B. EVALUATING TEST CELL PARAMETERS 

In the previous section of this chapter a relationship for 6, was obtained in terms 
of the fringe shift, AF , full width at half maximum, IT, and some test cell parameters, 
equation 4,23. Before determining the density profile, these parameters were evaluated. 

First consider the parameter 4 which represents the wavelength of the laser used in 
the interferometer. This experiment used a HE-NE laser for the interferometer light 
source Which has a wavelength of 633 nm in vacuum. This wavelength changed with 
each gas mixture used, however it was determined that the change, Inm, was not sig- 
nificant. Hence the value of / was fixed at 633nm throughout these calculations. 

Next, consider the parameter, p,, this represents the density of the gas im thewmen 
ence leg of the interferometer. It was assumed that the ambient conditions for each ex- 
periment Were the same, 300K and 600 torr. From the ideal gas law the ambient density 


was found to be: 


P 600/760 


a aes 3 . 19 3 
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The third parameter, K , is proportional to the molar refractivity of the gas mixture. 
(See equation 4.12). The following describes the procedure used to obtain the molar 
refractivity and subsequently the value of K for each gas mixture. 

Born and Wolf state that the molar refractivity of a gas mixture is equal to the 
weighted average of the molar refractivities of the constituents [Ref. 14: p. 89]. For the 


NH3-N2 experiment this can be written as 


Nya34dny3 t+ Nyd v2 


A Nv — (4.24) 
NH3—A2 : : 
Nyy3 tN a2 


where 
Sine, 1S the molar refractivity of the gas mixture 
ea, 1S the molar refractivity of NH3 
e A, is the molar refractivity of N2 
oe iS the numoer Of NHS molecul¢s in the test cell 


Se is tie numoer Of NZ miolecules in the test cell 


Therefore to obtain K for a particular gas mixture, the molar refractivity of NH3 
and N2 must be determined. The molar refractivity, equation 4.7, 1s a function of the 
index of refraction and the density. 

Consider first molecular nitrogen, N2. In order to determine its molar refractivity, 
A, the index of refraction at 633 nm was determined. This was accomplished by in- 
terpolating between the values given for 546.1 nm, 1.0002998, and 656.3 nm, 1.0002982, 
at 0° C and 760 torr [Ref. 15}. The value obtained, 1.00029854 was then corrected to 
300K and 600 torr using the Biot-Arago law [Ref. 16: p. 533}. 


eo P 
~ 1+ .00367¢ 760 


co 


where 
e n,is the index of refraction of a gas at temperature, t (degrees Centigrade) 
See is the index of refraction of a gas at 0’ C and 760 torr 
® tis the temperature of the gas in degrees centigrade 
See is tie ambient pressure in torr. 
Having performed the above corrections, the index of refraction for N2 at 300 K and 
600 torr was determined to be 1.0002 14440. 
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The tabulated value of the index of refraction for NH3 at the wavelength of sodium 
D light, 589.5 nm, was available from the Smithsonian Tables [Ref. 16: p. 533]. Insuffi- 
cient information prevented interpolating to obtain the value at 633 nm. Considering the 
difference to be quite small, the value at 589.5 nm, 1.000376, was used throughout these 
calculations. This value was then corrected for temperature and pressure using the same 
procedure outlined previously for N2. The index of refraction at 300K and 600 torr was 
determined to be 1.0002700. 

The next step in determining the molar refractivity was straightforward application 
of equation 4.7. This vielded a molar refractivity for N2 as 4.517 cm} per mole and for 
INS 5-088 cur per nigle. 

Having determined the molar refractivitv for N2 and NH3 equation 4.24 was used 
to determine the molar refractivity of the gas mixtures. These values are listed in Table 


2 along with the values of K , which were determined from the previously developed re- 





lationship 
3.4 
K= me (4.11) 


As will be shown later, K is the link between the Gaussian and Abel techniques. 
Having determined the numerical values of 4, p and K, the relationships for density, 
equations 4.16 and 4.23, can be reduced to the following: 


For uniform square channel 


col —— (4.26) 


For non-uniform square channel 


; _ _AF(O) 
a 1.0596P 


(4.27) 


K 
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Table 2 lists the value of @ and other mportant parameters for each gas mixture 


used in the experiment. 








@Geeeesl CASE 

Thus far in the development of the Gaussian analysis an algebraic expression has 
been obtained which allows one to determine the ratio of on-axis density to ambient 
densitv for a non-uniform channel, equation 4.27. A test case was generated to check 
the accuracy of this expression prior to applying it to the actual interferograms. 

The test case was a hypothetical cell filled with 600 torr N2 at 300K. Equation 4.21 
was used to develope a fringe pattern in which [ was equal to the square root of 2.8 and 


6, chosen to be 0.1. These conditions reduced equation 4.20 to the following: 
AF(Y) = 5.40 exp( —¥°) (4.28) 


Equation 4.28 was then solved for the value of Y which vielded one-half the maxi- 
mum fringe shift. Twice this value, 1.67 cm, corresponds to the full width at half maxi- 
mum. Evaluating equation 4.27 at this value of [, yielded a channel density of 0.10. Thus 


the equations developed provided self-consistent results. 


D. APPLICATION TO EXPERIMENTAL DATA 
Having determined that the equations developed for a non-uniform channel pro- 
vided self-consistent results, each fringe line was analvzed. Due to the asymmetry in the 
actual interferogram, (See Figure 11), three different values of the full width at half 
maximum were obtained. These different values correspond to 
e I, the value obtained if the right side were the same as the left side. 
e I, the value obtained if the left side were the same as the right side. 


e [; the value obtained when the left and right side values are averaged. 


A coniputer program was developed which determined the value of T and the value of 
AF and then computed the value of the channel density. The results of this program are 


listed in Tables 3 through 8. 


E. SUMMARY 

The characteristics of the reduced-density channel are determined by the channel 
density and the full width at half maximum at a particular location along the test cell. 
These characteristics have been determined from a method of Gaussian analysis devel- 
oped by Czuchlewski and are listed in Tables 3 through 8. While this method assumes 
the density profile is Gaussian, it provides a rapid means of determining an approximate 


density. Recalling equation 4.27 


: AF{0) 


ae 1.0590r (4.27) 


it is clear that upon measuring the peak fringe shift and the full width at half maximum 
a rapid determination of the channel density can be made. 
In the next chapter the Abel technique will be used to obtain the same reduced- 


density channel characteristics. 


oe) 
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V. ABEL ANALYSIS 


In the previous section a method for determining the ratio of on-axis density to 
ambient density and width of the reduced density channel was developed assuming the 
density profile had a Gaussian shape. In this section a method will be developed which 
makes no such assumption concerning the density profile, but requires the profile to be 
cylindrically symmetric. This technique was originallv developed by Neils Heinnik Abel 
in the early 1800's in his solution of a mechanics problem. In this section the Abel 
technique will be developed, applied to a test case and then applied to the interferometric 
data obtained during the reduced density channel experiment involving NH3-N2 gas 
nuxtures. The results of this analvsis are listed in Tables 9 through I4 and provide a 


determination of the density profile which is representative of the interferometric data. 


A. ABEL INTEGRAL TRANSFORMATION 

In his short life Neils Heinnk Abel made numerous contributions to the field of 
mathmatics and science. The significance of many of his contributions were not recog- 
nized until more than 100 years after his death in the late 1820's. One area of science 
Which has greatly benefitted from his work is the study of plasma spectroscopy. 

While attempting to solve a mechanics problem regarding the velocity of an object 
along an arc [Ref. 17: p. 97], Abel determined that ifa given function, H(a), can be re- 


lated to another function ‘¥ by 


| (5.1) 
0 (a— xy 


then ‘¥ can be obtained from 


x 
Y= ae it (5.2) 
9 (x-a) 


where 0 < uw < 1. The full developinent of this result is available in pages 97 through 


103 of Reference 17. 


As mentioned this integral transformation 1s widely used in plasma spectroscopy. In 
that apphcation the enussivity of a plasma source is determined from the mtensitv 
measurements. In the development of this emissivity-intensity relationship for use in 
plasma spectroscopy. Griem [Ref. 18: p. 176] found the intensitv of an optically thin 


plasma contained by a cylindrical column of radius R, to be 


Pe 


KY) =2 fv dX = io 


ae (333) 
: (r° a Ye 


Using Abel’s integral transformation the following determination of the emissivitv was 
obtained (See Figure 12) 





R 
Sl a ere) : 
é(r) = —; \ Ty (YS) ar (5.4) 
ali yee ee 
where qv the derivative of [(Y) with respect to Y. This same method is used 


throughout plasma physics in order to determine the emissivity of the source (See Refs. 
19 and 20 for examples). 

Typically the intensity is a measured quantity and not a differentiable function, 
hence the solution is accomplished bv numerical analvsis. Numerous methods of solving 
equation 5.4 numerically have been identified, some of these may be found in references 
21 through 25. 

Having identified the Abel integral transformation and one of its more widely used 
forms, the next step was to find an expression between the fringe data and the density 
which satisfied the relationship expressed in equations 5.1 and 5.2. Once identified, the 
Abel transformation was performed and resulted in an expression for the density as an 


integral function of the fringe shift. 


B. APPLYING THE ABEL TRANSFORMATION 
In the Gaussian development an integral expression between fringe shift and density 
was obtained for the case of the non-uniform channel, equation 4.18. Considering the 


radial symmetry of the gas, defined by 


Peay (24 
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Equation 4.18 was written as 





AF(Y) = 


Ion eo 2. : 
[_f80 ~ ax (3.6) 


Substituting for dX, its equivalent vielded 
Re 
2Kp, (d(r) — l]rdr 


AF(Y) = Sa ely? (S27) 


Performing the Abel Transformation on equation 5.7 yielded an expression for ratio of 


the channel density to ambient density 


R 

: ee ee ee 2) 5 

(ir) -1l= aL WY (y°—ry "dy (5.8) 
AE Ye, ee ; as 

aeve 15 the derivative of AF with respect to Y. Flence the ratio of the densities, 


dY 
O(r), can be determined by solving the right side of equation 5.8. Attempting to solve 


this analytically for any but the simplest cases ts extremely difficult. Consider as an ex- 
ample the Gaussian test case of the previous chapter. Substituting the derivative of the 


Gaussian fringe pattern into equation 5.8 produces 





R 
10.8 Y exp( — ¥7)\dY 
ek al (5.9) 
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Thus even for the situation where the fringe data is given as a differentiable function, 
numerical analysis is required to solve the right side of equation 5.8. 

Having determined that the Abel Integral transformation can be applied to the case 
of the reduced-density channels, the next step is to determine the method used to obtain 
a numerical solution to equation 5.8. 

As mentioned use of the Abel integral transformation in determining the emussivitv 
ofa plasma source is well documented. It is believed that this is one of the first uses of 
this method in determining the density profile from interferometric data. Solution of this 


expression, like most Abel integrals is best performed by numerical analvsis. In this case 


an analytic solution is not readily available. In the case of plasma spectroscopy the in- 
tensity is a measured quantity, and not a differentiable function. The next section de- 
scribes the method used to determine the reduced density channel characteristics using 


a numerical solution to the Abel integral. 


C. NUMERICAL SOLUTION 

A sketch of the test cell is provided in Figure 13. This slice of the test cell provides 
the foundation upon which the Abel technique was developed for obtaining the reduced 
density channel characteristics of the NH3-N2 gas mixtures. 

In the development, it was assumed that a radially symmetric density profile existed 
after the gas was energized bv the incident CO2 laser beam. Further it was assumed that 
the region could be divided into an infinite number of concentric regions in which the 
index of refraction was constant. To illustrate this Figure 13 1s divided into four regions. 
The index of refraction at all locations in region o is the same. Likewise the index of 
refraction for all locations in region y 1s the same and so on for regions y and v. This ts 
not to imply, however that the index of refraction in region a is the same as in regions 
y,4 or v. Having made these assumptions, an expression for the optical path difference 


along a chord, C,, was deternuned in terms of the notation in Figure 13. 
2 725 V; 2 meaner 2 es Z 
C= 2, 00? = 1) ee (5.10) 


Expanding this as a series of N regions yields the following 


N=1 


iz Y 
G=2) = ING2,=7)" ==) 
ji 


ey 


} (5.11) 


In order to obtain an expression for the index of refraction, the nght side of equation 


5.11 was manipulated to obtain 
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2 one 2 2a 
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In continuing with the development of the Abel technique, a relation between the 
index of refraction and density was required. This relationship was identified in the de- 


velopment of the Gaussian technique as equation 4.11. 
na—-l=Kp (4.11) 


Writing the previous equation in terms of density and substituting the equivalence of 


optical path difference, OPD = AF (Y) A, yielded the following 


N=} 


=i+] 
Pi ee 1 (5.13) 
2K(r2., — 17)” 





This equation can be solved by starting at the outer edge of the test cell, X =R, ,and 
proceeding towards the center, X=O. Having obtained an expression for p,, the next 
step was to determine the ratio of gas density in the channel to the ambient density. This 
ratio was obtained by scaling the density in the channel. Recall that the density along 
the axis was determined to be less than the density along the walls of the test cell. De- 


fining a scale comparable to that of the Gaussian technique, 6, was set equal to: 


(5.14) 


One can observe that from equation 5.14, the value of 6, is 1 at the wall of the cylinder, 
P-<p, , and decreases at positions close to the center of the channel. 

Having determined a relationship for the density in the channel, a computer pro- 
gram was developed which extracted the density profile of the channel from the fringe 
data. This method was then checked using the Gaussian test case developed for the 
Gaussian technique. Application of this test case to the Abel technique is described in 


the next section. 


D. TEST CASE 

A method of obtaining the density profile of the test cell using the Abel technique 
was developed in the previous section. This method was checked using the fringe pattern 
described by equation 4.28 and found to provide reults of 6, which agreed with test re- 
Suitssto Within two decimal places; the results of T agreed to within 0.07 cm. Table [5 


Peete Abel and test results for 6, and [ for each gas mixture used. As an example of 


the resultant density profile, the Abel determined prolile for 600 torr N2 obtained from 
the fringe pattern described by equation 4.28 is plotted in Figure I4. The values used 
to generate the plot are provided in Table 16. Having determined the Abel method 
produced results which agreed with the test cases, the density profile for each fringe line 


digitized was determined. 


E. APPLICATION TO EXPERIMENTAL DATA 

After determining that the Abel technique produced satisfactory results, each fringe 
line was analvzed to determine the density profile. In doing so it became clear that the 
slanted interferogram, characteristic of the NH3-N2 gas mixtures, presented a problem 
that was not encountered with the test case. 

In the test case, the fringe lines were divided into two parts and the density profile 
of each part matched at the origin. In the case of the actual interferograms the difference 
in fringe patterns on either side of the division line precluded a match at the origin. This 
difference is consistent with the development of the Abel technique. however it violates 
the assumption that the density profiles were cylindrically symmetric. At this point the 
assumption was made that each part of the fringe line represented half of a symmetric 
fringe line. The value of 6, to be used in the analysis of the results was an average of 
the values obtained for both parts of the fringe line. The same procedure was used for 
deternuning the value of [. The results of Abel analvsis on each fringe line are listed in 
Tables 9 through 14. The values of 6, are provided for both sides of the division line as 
Well as the average between the two sides. These three values can then be thought as 
defining the minimum, maximum and most likely value of the on-axis density and 


channel width. 


F. SUMMARY 

The Abel technique provided a means to determine the density profile of a gas 
mixture heated by a CO2 laser beam. From this density profile a determination of the 
reduced density channel characteristics was possible. In the next chapter results of the 
the Gaussian and Abel techniques will be used to determine relationships between the 


channel characteristics and the variables of the experiment, energy and NI13 pressure. 


VI. RESULTS 


Interferograms made shortly after a CO2 laser energized a test cell containing vari- 
ous mixtures of N2 and NH35 have been analvzed using two different methods. These 
methods, Gaussian and Abel, have been described in detail tn the preceding chapters. In 
this section a comparison of the two methods and some general observations will be 


made concerning channel formation. 


A. GAUSSIAN-ABEL COMPARISON 

A summary of the experimental results is provided in Tables 17-22. Included in 
these tables is a comparison between the Gaussian and Abel determined values. The 
closest agreement between the two methods can be found tn Table 19. In that Table the 
results of the experiment for a test cell contaming a gas mixture of 15 torr NII3 are 
listed. At 112 cm along the test cell and for a laser energy of 84 joules, the Gaussian and 
Abel determined values of 6, differ by 1.7 percent. Under the same conditions the values 
determined for I’ differ by 3.6 percent. This fringe pattern and the Gaussian curve ob- 
tained by using values of AF and I obtained from the fringe data are plotted in Figure 
15. As expected the two methods vield similar results when the fringe pattern has a 
Gaussian shape. 

In some cases, however, the fringe data are not Gaussian shaped and the analvtic 
analysis is not expected to be very accurate. Consider the results of 5 torr NH3 listed 
m Table 17. At 51 cm and for an energy of 31 joules, the values of 6, and [’ obtained 
by the Gaussian and Abel methods differ appreciably. For 6, the two analysis methods 
Mieleresults which differ by 21.1 percent; the values for I differ by nearly 331 percent! 
This fringe line and the corresponding Gaussian curve are plotted in Figure 16. It ts 
apparent from the differences between the curves that the fringe pattern is not Gaussian 
shaped. Recalling that the Gaussian method uses the maximum [fringe shift and the [ull 
width at half maximum to determine the value of 6, and that the Abel method uses all 
the points digitized along the [Tringe line, it ts easily understood why the two methods 
yield different results for non-Gaussian shaped fringe patterns. 

In general it can be stated that the smaller the difference between the Gaussian and 
Abel determined values of 6, and I the more nearly the fringe line approximates a 
Gaussian shape. A review of Tables 17-22 reveals that the Gaussian determined PF and 


fem wecl determined 1 differ b¥ more than 20 percent in all but 8 of the fringe lines. 


The corresponding values of 0, in these instances differed bv less than 3.3 percent: 
Hence less than 20 percent of the experimental data was Gaussian shaped. For this 
reason analysis of the reduced-density channel formation will be based on the results 
obtained from Abel method. 

The Gaussian method was developed to provide a quick assessment of the channel 
formation. From these experimental results it can be stated that for fringe patterns of 
Gaussian shape, it vielded results very similar to those obtained from the Abcl method. 
For non-Gaussian shaped fringe patterns, it can be stated that it vields values of 6, 
within 31 percent of the Abel method. As for the values of [ obtained, the Gaussian 
method provided results which were always greater than or equal to the Abel results. 
Hence, for this experiment the Gaussian method identifies the largest channel width 


possible. 


B. ACCURACY OF RESULTS 

In order to provide a determination of the accuracy of the entire process from pho- 
tographic enlargement through the determination of 6,, a series of error estimates Were 
made. In these calculations it was assumed that the largest probable error in digitizing 
a fringe line would be 10 percent of the actual value. Using the fringe pattern described 
by equation 4.28 two new fringe patterns were generated each differing from the original 
by 10 percent. One of the fringe patterns generated consistently differed by 10 percent 
more than the original pattern; the second pattern generated consistently differed by 10 
percent less than the original pattern. The results of the Abel deternuned values of 6, at 
each gas nuxture are plotted in Figure 17. As can be observed a difference of + 10 
percent in the fringe pattern results in a + 0.09 differnce in 6,. Since the assumption of 
{OQ percent error incorporates the error introduced in the photographic enlargement 
process and covers crratic digitizing errors, there 1s confidence that the values of 0, are 


accurate to within + 0.09. 


C. ANALYSIS OF RESULTS 

The Abel method is more dependent on the shape of the entire fringe line than the 
Gaussian method. To gain insight into the meaning of the results listed in Tables 17-22, 
several plots of the Abel determined values of 6, and F were made. The results could 


be plotted in numerous Ways, for example 0, as a function of NI13 concentrations 
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energy, distance from the entrance of the cell, etc. The methods chosen to represent the 


results are 6, as a function of laser energy, Figures 18-22, and 0, as function of distance 


from the test cell entrance, Figures 23 and 24. These two methods will provide a means 


to compare the results to the results of the SF6 reduced-density channel! experiment. 


1. Channel Density 


The following observations were made from the plots of 6, as a function of en- 


ergv and with NH3 concentration and distance from cell entrance constant parameters 
(Figures 18 through 22). 


Pam ecm itommtne entrance oO. the test cell{ See Figures 18 and 19); 


The plots for 5 torr and 200 torr NH3 indicate that there is a slight decrease in 0, 
as energy is increased. One could also fit the data with a constant value for 6, of 
0.73 + 0.05 for 5 torr and 0.83 + 0.03 for 200 torr. Hence an equallv valid inter- 
pretation of the results is that there is no apparent energy dependence on the 
channel depth (Figures 18 and 19). 


For 15 torr (See Figure 18) a pronounced dependence of 6, on the energy is ob- 
served, channel density decreases with increasing energy. At low energies 0, is ob- 
served at values around 0.60, while at higher energies values around 0.30 are 
observed. 


For 50 torr the channel density decreases with increasing energy through the range 
of 0 to 43 joules. At 43 joules there appears to be a minimum value and then the 
channel density begins to increase. This minimum will later be compared to the 
minimum predicted by the computer model for the SF6 experiment. 


Eom eciniitoim tie entralice to the test cell ( See Figures 20 and 21): 


For 5 torr there appears to be little effect on 6, when energy is increased. The re- 
sults could be fitted with a constant 6, of 0.73 + 0.07 (Figure 20). 


For 15 torr a distinguishable dependence of 6, 1s observed, channel density de- 
creases as energy increases (Figure 20). 


For 50 torr one can observe remnants of the minimum identified at 51 cm (Figure 
Din). 


ee cm ftom the entrance to tlre test cell (Sce Figure 22): 


For 5 and 15 torr NEH3 a slight decrease in channel density is observed with in- 
creasing energy. 


For 50 and 200 torr NH3 all the laser energy is absorbed at the beginning of the 
cell and thus, there are no channels at 221 cm. 


The results obtained for 0, are consistent and viable. The following identify traits 
supporting thus 


Woes Ol ©) versus encereyv (Fieures 18 throual 22) vield smooth curves: there are no 
discontinuities. 
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e The values obtained for 6, under similar conditions vield similar results. An exarn- 
ple is the results obtained for 15 torr at 51 cm and S4 joules (Figure 18). The values 
of 6, are close and are within the statistical error bar. 


To provide further comparisons, representative data points were selected and 
plotted as a function of position along the test cell in Figures 23 and 24. The energies 
were selected to provide a means of comparison between different concentrations of 


NH3. Among the more prominent features of each plot in Figure 25 are 


e As the distance from the entrance of the test cell is increased the value of 6, in- 
creases. This finding agrees intuitively with the idea that the laser beam is contin- 
uously absorbed as it travels further from the entrance. Hence, there is less energy 
available at the end of the channel than at the beginning and we expect weaker 
channels at the end of the cell. The situation for 5 torr, 45 joules and at distances 
less than 130 cm from the entrance of the cell mav be more complicated. 


e As the concentration of NH3 increases, the length of the channel decreases. This 
is evidenced bv the fact that no channels exist in 200 torr at 112 cm and none exist 
at 221 cm for 50 torr. It also agrees intuitively with the idea that as the number 
of absorbers per unit volume increases, the distance a beam travels before dissi- 
pating decreases. 


e As the concentration of NH3 increases from 5 to 15 torr, channel density 
decreases. This is true for all energies plotted and is consistent with the findings of 
the N2-SF6 experiment. In short, the higher the concentration, the deeper tlie 
density channel (Figure 26). 


e For a laser energy of 45 joules the channel produced in 5 torr NH3 is relatively 
constant over the range of 50 to 221 cm with a value of 6, of 0.79 © OOS, igus 
also true for 75 joules, for which 6, remains 0.68 + 0.04. This 1s similar to the re- 
sults obtained for the lowest concentration of SF6 [Ref. 8: p. 27], for which rela- 
tively long shallow channels were observed. 

2. Channel Width 
The channel widths for each concentration at selected energies are plotted in 
Figure 24, The selected energies are the same used in Figure 23 and allow a a comparison 
between values at the same energy. Unlike the plots of 6, versus distances, very little 
consistent information can be garnered from Figure 24. One identifiable characteristic 
can be observed by considering the 5 and 15 torr concentrations at 51 cm anduiiZseme 
In these instances the channel width increases as energy increases. The opposite appears 


true for 200 torr NH3, D decreases as energy increases. 


D. COMPARISON WITH SF6 

At present there is no quantitative method for comparing the characteristics of the 
reduced-density channel obtained in the NH3-N2 experiment. For that reason the results 
will be compared with the results of the N2-SF6 reduced-densitv channel experiment 
acies: p. 27}. 

The results of the N2-SF6 and the N2-NH3 experiments were very similar. The 
major differences being that NH3 was substituted for SF6 in the test cell as the 10.6 
um radiation absorber and the measurements were taken at different distances. A com- 
puter model, CHANL, developed by Naval Research Laboratory scientists was used to 
predict the results of the N2-SF6 experiment [Ref. 8: p. 8]. A signficant product of that 
model is plotted in Figure 25, This figure shows the predicted channel density for a laser 
energy of 50 joules as a function of distance from the entrance of the test cell. The 
\2-SF6 expernunental results, shown in Figure 26, correlate well with the results of the 
model. A trend of increasing absorber concentration resulting in shorter, deeper channels 
1S apparent in both plots. When comparing Figures 25 and 26 with the N2-NH3 results 
shown in Figure 23, the same pattern clearly exists for 5, 15 and 50 torr NH3. Less 
obvious 1s the correlation with 200 torr NH3. While only two data points are available, 
one can assume that the laser energy has been almost totally absorbed before it reached 
51cm. This assumption is supported by Figures 26 and 27. 

Figure 27 is another product of the CHANL code. This plot shows the effect on 6, 
at 0.4 torr SF6 as energy is increased. Significant is the minimum 6, which occurs at 50 
joules and the sharp change in slope at 165 cm. Evidence of these changes in slope are 
evident in both the N2-SF6 (Figures 28 and 29) and N2-NH3 (Figures 19 and 21) ex- 
perimental results. 

Based on the above observations, it can be concluded that the N2-NH3 results ex- 
hibit the same qualitative behavior as the N2-SF6 results and the CHANL code pred- 


ictions. 


E. CONCLUSIONS 

The CO2 laser beam created reduced-density channels in test cells containing a gas 
mixture of N2 and NH3. These channels were analyzed and found to be dependent on 
many parameters including the concentration of NH3 and the laser energy. As expected 
femassound that the depth of the channel decreased as the the beam travelled along the 
test cell. It was also determined that for a fixed laser energy. the depth of the channel 


and its length varied inversely with the concentration of NII3. Vhat is. shorter, and 
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deeper channels were observed at higher concentrations of NH3; long, shallow channels 
were observed at lower concentrations. [his qualitatively agrees with the experimental 
results and computer predictions for the N2-SF6 experiment [Ref. 8: p. 27]. 

It can also be concluded that the Gaussian method of analvsis provides a quick 
means of determining the characteristics of a reduced-density channel. The accuracy of 
the method is dependent upon the shape of the fringe pattern. For Gaussian shaped 
patterns the results are highly accurate. For non-Gaussian shapes values for 0, differed 
from the more accurate Abel method by as much as 31 percent. The channel width de- 
termined differed bv as much as 330 percent. Despite these large differences the Gaussian 
method provided a quick means of determining 6, to within .15 and finding the maxi- 
mum value of channel width. 

To summarize, the analysis of reduced-density channels in gas mixtures of N2 and 
NH3 support the N2 and SF6 experimental results. The channels obtained in both ex- 
periments should be suitable for studying the propagation of intense electron beams 
through the atmosphere, and thus provides a foundation for determining the role these 


beams will have in the Strategic Defense Initiative. 
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Table 1. COORDINATES OF ACTUAL DATA AND ENLARGED DATA 
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Table 3. GAUSSIAN RESULTS FOR 5 TORR NH3/395 TORR N2 
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Table 4. GAUSSIAN RESULTS FOR 15 TORR NH3/385 TORR N2 
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Table 5. GAUSSIAN RESULTS FOR 15 TORK Ni3/333 7ORieee 
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Table 6. GAUSSIAN RESULTS FOR 30 TORR NH3/550 TORR N2 





























FRINGE Dts ENE KG 


(joules) 


GAUSSIAN ANAL ISIS 
(ONeeeSIS DENSITY) 






a eS a ce 
jTest_ [i NANA aa 
ps77ao_[o S13 fas S352 83 
smas |s_|st ie —~id asd ao—ids d 
a Psa: 76 SSC*YSCT <a 90) amas 72 738 
revsps_ [8 | fe —~dt to ~<a) to 
sao [ost fe SIS? ST SB 
a 
io 2 eS Saaee ee Se 3 

rssop7 [7s —~*d8 
ac 





tr 
cn 
~~] — 


a 
cn 






tt 

ae Cn un 
oa Carp to 
i. j— 


> 
a, 
> 


4 


Table 7. GAUSSIAN RESULTS FOR 30 TORR NH3/550 TORR N2 
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Table 8. GAUSSIAN RESULTS FOR 200 TORR NH3/400 TORR N2 
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Table 9. ABEL RESULTS FOR 5 TORR NH3/595 TORR N2 
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Table 10. ABEL RESULTS FOR 15 TORR NH3/385 TORR N2 
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Table 11. ABEL RESULTS FOR 15 TORR NH3/385 TORR N2 
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Table 12. ABEL RESULTS FOR 50 TORR NH3/550 TORR N2 
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Table 13. ABEL RESULTS FOR 50 TORR NH3/550 TORR N2 
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Table 14. ABEL RESULTS FOR 200 TORR NH3/400 TORR N2 
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Table 16. DENSITY PROFILE OF 600 TORR N2 (ABEL) 
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Table 17. SUMMARY OF RESULTS BOR STO Khai 
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Table 18. SUMMARY OF RESULTS FOR 15 TORR NH3 
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Table 19. SUMMARY OF RESULTS FOR 15 TORR NH3 
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Table 20. SUMMARY OF RESULTS FOR 50 TORR NH3 
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Table 21. SUMMARY OF RESULTS FOR 50 TORR NH3 
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Table 22. SUMMARY OF RESULTS FOR 200 TORR NH3 
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Schematic of Reduced-Density Channel Experiment 


mieure 1. 
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Figure 2. CQO2 Laser Beam Pulse Characteristics 
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Fisure 4. Typical {nterferogram after Energizing CO2 Taser 
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Figure 5. Sketch of Fringe Lines with Key Characteristics Annotated 
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Figure 6. Pattern Used to Test the Enlargement Process 
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Sketch of Interferogram Annotated to Resolve Geometry 


Figure 3. 
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Sketch of Interferogram from Square Channel 


Figure 9. 
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Sketch of Interferogram from Cylindrical Gaussian Channel 


Figure 10. 
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Figure 11. Sketch Showing [ of Asymmetric Fringe Line 
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Figure 12. Schematic View of a Cylindrical Plasma Column 





Figure 13. Cylindrically Symmetric Test Cell Containing NH3-N2 Gas Mixture 
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Abel Determined Density Profile of Gaussian Fringe Pattern 


Figure I4. 


66 


(9) FONVLSIG WIGVHY 
a G'} | G 0 G | Git 





WwoZLL=Z fps = "Q "HN wOlSL :vLVG SONIHA 


ee oes islet 


Comparison of Fringe Data to Gaussian Curve; Best Case. 


Figure 15. 
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Figure 16. 
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Comparison of Fringe Data to Gaussian Curve; Worst Case. 
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Figure 17. 
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Plot of On-Axis Density (Abel) vs NH3 Concentration for Test Case and 
Test Case + 10 Percent. 


69 


Z =o | China Oleh ime 


1.0 





20 30 40 50 60 70 80 
ENERGY (Joules) 


Z2= 51 Cilio (Chenin 





20 30 40 50 60 70 ~—«-80 
ENERGY (Joules) 


Figure 18. Variation of Channel Density with Laser Energy at 51 cm from Entrance 
of Test Cell for 5 torr NH3 (top), 15 torr NH3 (bottom). 
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Figure 19. Variation of Channel Density with Laser Energy at 51 cm from Entrance 


of Test Cell for 50 torr NH3 (bottom), 200 torr (top). 
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Figure 20. Variation of Channel Density with Laser Energy at 112 cm from En- 


trance of Test Cell for > torr NE3 (top). 13 tor iio (uot om: 
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Figure 21. Variation of Channel Density with Laser Energy at 112 cm from En- 


trance of Test Cell for 30 torr NH3. 
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Figure 22. Variation of Channel Density with Laser Energy at 221 em from En- 


trance of Test Cell for 5 torr NHS (top), [3 ton i Doon: 
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Variation of Channel Density with Distance from Test Cell Entrance for 


mremre 23. 


0 and 200 torr NH3. 
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Variation of Channel Width with Distance from the Test Cell Entrance 


Figure 24. 
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Predicted Channel Density for 30 joule Laser Energy for Concentrations 
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Figure 26. Variation of Channel Density with Distance from the Test Cell Entrance 
for 0.2, 0.4, 03s and 220 torr Sire 
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Figure 27, Predicted Variation of Channel Density with Laser Energy for 0.4 torr 
SF6 at 50, 165 and 220 cm from the Test Cell Entrance. 
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Figure 28. Variation of Channel Density with Laser Energy for 0.4 torr SFO at 17 
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cm from the Test Cell Entrance. 
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Figure 29. Variation of Channel Density with Laser Energy for 0.4 torr SF6 at 165 


cm from the Test Cell Entrance. 
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